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Abstract 

We apply the methods of Field Theory to study the turbulent regimes of statistical 
systems. First we show how one can find their probability densities. For the case of the 
theory of wave turbulence with four-wave interaction we calculate them explicitly and 
study their properties. Using those densities we show how one can in principle calculate 
any correlation function in this theory by means of direct perturbative expansion in powers 
of the interaction. Then we give the general form of the corrections to the kinetic equation 
and develop an appropriate diagrammatic technique. This technique, while resembling 
that of <^> 4 theory, has many new distinctive features. The role of the e = d — 4 parameter 
of (p 4 theory is played here by the parameter k = (3 + d — a — 7 where f3 is the dimension 
of the interaction, d is the space dimension, a is the dimension of the energy spectrum 
and 7 is the "classical" wave density dimension. If k > then the Kolmogorov index is 
exact, and if k < then we expect it to be modified by the interaction. For k a small 
negative number, a < 1 and a special form of the interaction we compute this modification 
explicitly with the additional assumption of the irrelevance of the IR divergencies which 
still needs to be verified. 
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1. Introduction 

The theory of wave turbulence studies the stationary states of the statistical classical 
(not quantum) system consisting of waves with a small interaction, (on the theory of wave 
turbulence see ref. |jj and references therein) . Its Hamiltonian can be written down in the 
following form 

P2P3P4. a pi a p2 a P3 a P4 (1-1) 

P P1P2P3P4 

It is just a collection of waves with the energy spectrum uj p and the four wave inter- 
action A Pl p 2 p 3 p 4 with the evident properties A Pl p 2 p 3 p 4 = Ap 2 p 1 p 3 p 4 = Ap 3 p 4 p lP2 . Let us note 
that this Hamiltonian conserves the total wave number 

N = j24 a p- ( L2 ) 

p 

The simplest stationary state of this system is a thermodynamic equilibrium and its 
probability density is given by the well-known Gibbs distribution exp(— H ~tp N ). The basic 
property of thermodynamic equilibrium is that the detailed balance principle is satisfied. 
There are as many waves going from the wave- number pi to p2 as there are ones going 
back. However there are other stationary states where that principle is not satisfied, or in 
other words, while the total number of waves coming to the given wave number p is zero, 
there is a flux of waves through the system. Which state will be chosen by the system 
depends on the external conditions. If it interacts with a heat bath satisfying the detailed 
balance principle, it will soon settle into the thermodynamic equilibrium. If, on the other 
hand, the heat bath is so special that it injects waves to the system at one wave number 
and removes them at a different one, then the system will necessarily choose one of those 
extra states. A simple example of the latter case is the waves on the surface of water which 
are injected by, for example, the ship, at wave lengths of the order of ship length, and are 
dissipated at much smaller lengths by viscosity. 

The theory of turbulence concerns itself with studying the "flux" states as the thermo- 
dynamic equilibrium has already been studied in great details by Statistical Physics. The 
standard hydrodynamics turbulence is also the example of the flux states since here we 
have a very complicated motion of liquid with a stationary probability distribution which 
is characterized by the flux of energy or other conserved quantities through various scales 
of vertex motion. 
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While the existence of those extra states was known for a long time and they were 
analyzed with the aid of kinetic equations (ref. H), the explicit form of the probability 
distribution was never found. But let us remember that the Liouville theorem of the 
statistical mechanics tells us that the probability distribution of any system should be 
an additive integral of motion. The next step is usually to assume the system doesn't 
have any additive integrals of motion different from the energy and the particle (or wave) 
number. That's how Gibbs distribution is obtained. But in fact any dynamical system 
has as many integrals of motion as the number of its degrees of freedom. Really, any 
initial conditions expressed in terms of the changing variables give us a formal integral of 
motion. The important role such integrals play in the description of the flux-states of the 
statistical systems was first noticed by A. Polyakov (ref. @). Essentially such probability 
distributions mean that we prepared first a statistical ensemble with the probability p(q,p) 
where q and p are the initial coordinates and momenta respectively. Then we follow its 
time development. A correlation function of some quantity X(q,p) at the moment t = 
is given by 



To find this function at a later time t we need to use the solutions of the equations of motion, 
p(t) and q(t), to express the initial conditions in terms of the values of the coordinates and 
the momenta at a later time t, q(q(t),p(t),t) and p(q(t) , p(t) , t) . The correlation function 
will then be given by 



where p' is the same probability distribution expressed in terms of the current coordinates 
and momenta q(t),p(t). The Liouville theorem of classical mechanics dpdq = dp(t)dq(t) 
allowed us to make that change of variables. 

Constructed in this way, p' is not a stationary probability distribution. But we expect 
it to have some limit as t — > oo otherwise the system we study does not have stationary 
states. But even if such a limit exists, it does not automatically imply all the correlation 
functions will also have some well defined limit as t — ► oo. This phenomenon will be dis- 
cussed later in this paper. Now it will be enough to say that the additional conditions p' 
should satisfy to give a stationary state for the system is called the kinetic equation. To 
uniquely find the probability distribution one must solve the kinetic equation. Thermody- 
namic equilibrium is often just one of many solutions. 





3 



The main problem of the theory of turbulence is to solve the kinetic equation and 

v v 



to find the correlation functions, especially the two point correlation function < a) a p > 



which gives the spectrum of the theory. 



In this paper we are going to study the kinetic equation for the system ( |1.1| ) and 
develop a diagrammatic technique to represent the infinite series in powers of A we will have 
to deal with. This technique while somewhat similar to that of <^> 4 theory will have many 
features not encountered in Field Theory before. It may be even viewed as a generalization 
of if technique since in the limit of thermodynamic equilibrium it reduces to the standard 
<£> 4 theory. 

Nevertheless we are going to find some specific cases, namely when the interaction 
coefficient A is just a product of the external momenta, when one can extract the infinite 
series of most divergent diagrams, just like in <£> 4 theory in four dimensions, with the 
condition of the dimensionlessness of the interaction coefficient, of course. With some 
additional assumptions, in particular neglecting the IR divergencies, we shall make use of 
it to calculate the correction to the "classical" spectrum of the theory. 

Technically there is no need to know the probability distribution to compute all the 
correlation functions of the theory ( |1 . 1| ) as discussed later in the paper. But we hope 
that the study of the probability distribution for will allow us to understand better 
the behavior of the correlation function of the hydrodynamic turbulence where no reliable 
methods of finding the correlation functions exist. We shall discuss that in the end of the 
paper. 

The method which was widely used to compute the correlation functions thus far is 
called the Wyld's technique (see refs. @, Q). It is essentially a classical nonequilibrium 
diagrammatic technique, and it deals with the development of the correlation functions in 
time unlike the technique proposed in this paper which deals with the equal time correla- 
tion functions, as in statistical mechanics. Recently a paper || was published where the 
authors claimed they could in principle compute the kinetic equation up to any order of 
the interaction using the Wyld's technique. We think that while the kinetic equation ob- 
tained in this paper should be equivalent with that of || , our method is easier for practical 
calculation. 



2. The Probability Distribution 

According to the program outlined in the Introduction, first of all we need to study 
all possible integrals of motion for the system ( |1 . 1| ) . Let us see how one can construct 
those integrals. We shall start by assigning the variables a p and a| the initial values 
a p (t = 0) = dp, a p (t = 0) = a p °. Then by solving the equations of motion one can find 
a(t) and a^(t) as functions of a and and time. By inverting those functions, one can 
find a and as functions of time, a(t) and a^(t). But by definition a and do not 
depend on time. So they are the integrals of motion. If we want to construct something 
out of them which can play the role of the density of waves, we should consider the linear 
combination 

p 

where f p are some arbitrary coefficients. After being expressed in terms of a(t), a} (t) and 
t it becomes a valid integral of motion. Its explicit time dependence can be eliminated by 
passing to the limit t — ► oo. 

In order to find the explicit form for F we need to solve the equations of motion. The 
Hamiltonian (|1 . 1| ) allows us to solve them perturbatively and we can obtain F in terms of 
a series in powers of A. We can even avoid solving the equations of motion if we use the 
following procedure. One should look for F in terms of a power series 

F = ^ ] f p a p a p + ^ ^ PlP2P3P4 a jp 1 a p? a P3 a P4~\~ 

p (2.2) 

+ ^PlP2P3P4P5P6 a pi a p 2 a p3 a P4 a P5 a P6 + • • • • 

Here A and O are some still unknown functions. We impose the condition on F that it is 
an integral of motion, or {HF} = 0, { and } being the Poisson brackets. It allows us to 
find those functions to get 

A _ fpi ~l~ fp2 ~ /p3 ~ /p4 \ (r) o\ 

iV PlP2P3P4 — , , , , . , , , , • A PlP2P3P4 

u Pl + u P2 — u P3 — uj Pa — ie 

Q _ A ST~^ (-^P7PlP5P6^P2P3P4P7 ~ ^P7PlP5P6 ^P2P3P4P7 ) Icy A\ 

1 2 PlP2P3P4P5P6 / j i i o • ' \ Z -^J 

^ Vm + uj P2 + uj P3 - u PA - uj P5 - u P6 - 2it 

We can in principle find recursively all the terms in the series ( |2.2| ), one after another. 

A few words must be said about e's which appear in the denominators. Technically, 
when we compute Poisson brackets no e's appear. But we must introduce them to avoid 
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the poles in ( P-3| ) and ( j2.4j) in order for the sums in ( |2.2| ) to make sense. This procedure is 
legitimate only if we deal with the infinite number of degrees of freedom, that is the sums 
in ( |2.2| ) are really the integrals. If we had dealt with finite dimensional system, we would 
have obtained unavoidable zeros in the denominators of (|2.2| ) and the integrals of motion 
we discuss would not have existed. This is a very important observation. 

Another way to see why e must appear is to use the original procedure of finding F. 
One would discover that the function F depends on the expressions like (1 — exp(— uot))/u 
and the only way to pass to the limit t —>■ oo is to assume that u has an infinitesimal 
imaginary part. The last procedure ensures the signs of e's are correct. 

It may seem somewhat strange that we paid attention to writing He in the denominator 
of O as l/(x — ie) = 1/x + in5(x) regardless of the magnitude of e. However, later we shall 
see the cases when the magnitude of e becomes important, in particular when we deal with 
the expressions like e8(x)/(x — ie). To check the magnitude of e of the denominator of fi, 
one should follow the finite time procedure described in the previous paragraph carefully. 
One would discover that the denominator of O came from summing up the denominators 
of two A's included in the definition of O. If each of these A's has — ie in its denominator, 
then there should be — 2ie in the denominator of O. This point will be clarified later in 
the paper. 

The parameter e plays the role of a regulator in the theory to be discussed and in this 
sense is analogous to the cutoff used in Quantum Field Theory^. The function F by itself 
is only an approximate integral of motion as long as e is finite and we must pass to the 
limit e — > in all our calculations. 

The construction of the integrals of motion ( [2.2| ) represents a certain paradox. The 
system given by ( |1 . 1| ) is not necessarily integrable. Moreover, we are not interested in 
integrable systems. If the system were integrable, we could always make a change of 
variables to bring it to the form of noninteracting oscillators where no fluxes of energy 
could exist. Nevertheless we have just constructed the complete set of the integrals of 
motion for ( p. . 1|) . Something must be wrong with those integrals, otherwise the system 
( |TTT| ) would be integrable. 

It is those poles in the denominators of (|2.3| ) and (|2.4| ) which are responsible for those 
integrals not to give rise to any integrability. They actually make the integrals ( |2.2|) to 
be multi-valued functions which cannot be used to construct those tori the trajectories of 

' It is a time cutoff. The momentum cutoff will also have to be imposed if needed. 
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integrable systems lie on. But ( |2.2| ) can perfectly be used to play the role of probability 
distribution if one also introduces e's to regularize them. The finite time computation, as 
was discussed above, validates the introduction of e's in the form as in (|2.3|) and (|2] 



Now our statement is all the statistical states of the theory ( |1 . 1| ) can be described by 
the probability density given by 

exp(-F). (2.5) 

Even a thermodynamic equilibrium is a particular case of fl2.5|) , with f p = u p or f p = 1. 
All the correlation functions of the theory can be found using 

<X>=^ J H da p H da\X exp(-F) (2.6) 
p p 

where Z is a partition function 

Z= / Y[da p JJda+exp(-F) (2.7) 

p p 

In a complete analogy with Quantum Field Theory, the role of the partition function in the 
denominator amounts just to removing "vacuum diagrams" , or contractions totally inside 
the interaction, so we will usually omit it in our formulae while keeping that in mind. 

We are going now to demonstrate how one can derive a well-known kinetic equation 
(in the first nonvanishing order of perturbation theory) using our probability distribution. 
The kinetic equation is just a statement that the derivative with respect to the time of the 
correlation function < a p a p > is zero. In other words, 

d(a\.a v ) r . , 

at . (2.8) 

= ~~ * / j ^PlP2P3P4 ^ a pi a p 2 a P3 a P4 ^ {fippi ^PP2 _ ^PP3 ~ ^PP4 ) = ^ 
P1P2P3P4 



Since the function F is real, one can easily see from ( |2.8p that only the imaginary part 
of the correlation function < a pi a p2 a P3 a P4 > will contribute to the kinetic equation which 
can actually written down in the form 

-4i A pp 2 P3P4 Im < « P a P2 a P3 a P4 >= (2-9) 

P2P3P4 
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Let us find the four point correiation function perturbatively. Its zero order value is 
equal to 



— / J^da P Y\da p a pi a p2 a P3 a P4 exp( 

^ ] fp a p a p) 



p 

~~f 7 (^PlP3^P2P4 ^PlPA^P2Ps) 
JpiJl 



P2 

This expression is real so it will not contribute to the equation ( |2.9| ). Now the first order 
contribution to that function is obtained by expanding the exponent in powers of A and is 
given by 

_A\ fp3 + fp4 ~ fpi ~ fp2 1 , p (O m\\ 

^ A PXP2P3PA , _ _ _• | f f f l Z - iU J 

Lfp 3 T LUp 4 U/ pi 0Up 2 ^ J p\J P2J P3J P4. 

where R is a real expression coming from the contractions inside the interaction which is 
not interesting because it is real and also because it is zero in case if the interaction A 
conserves the momentum. Imaginary part of (|2.10| ) is 

4Z7T A pi p 2 p 3 p 4 {fp 1 + f P2 ~ fp 3 ~~ fp4,)fi(w Pl + uj P2 — uj P3 — iv P4 ) — - - - 

J Pi JP2JP3 J Pa 

and after substituting it to ( |2.9| ) we arrive at the kinetic equation 

ie - e x iv2P.pA"v - "P2 - - ^j /p+ / p r/ p r /p4 = °- ^ 

P2P3P4 JpJp2J P3 JpA 

One can easily recognize the standard kinetic equation of ref. if one substitutes 
with the same accuracy < a p a p >= 1/ f p . 



3. The Anomaly 

As it has already been mentioned in the Introduction, it is highly nontrivial that the 
kinetic equation we obtained is not identically zero. How can a stationary probability 
distribution give rise to a correlation function depending on time? More precisely we can 
use the following argument. For any quantity X 

< ^ >=< {HX} >= J {HX}exp(-F) = J {HF}Xexp(-F) = 

where we took the integral by parts and used that F is the integral of motion, that is 
{HF} = 0. But it contradicts our computations. The resolution of this paradox lies in 
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the fact that F is the integral of motion only in the limit e — ► 0. If we compute {HF} 
we shall discover that it is of the order of e. But then, when computing the integral 
J{HF}X exp(— F), we shall find that it in itself is of the order of 1/e. Those epsilons 
cancel and we can safely pass to the limit e — > to get the finite answer for the kinetic 
equation. The effect is completely analogous to the anomaly of Quantum Field Theory. 
By imposing the cutoff e on the theory we violate the conservation of F and it is never 
fully restored after removing the cutoff ^. 

It is instructive to rederive the kinetic equation in this fashion. The Poisson bracket 
{HF} can be recursively computed to give 

{HF} = —eAp 1 p 2 p 3 p 4 ap l ap 2 ap 3 a P4 — 2eflp 1 p 2 p 3 p 4 p 5 p 6 ap l ap 2 ap 3 ap 4 ap 5 a P6 — . . . . (3.1) 

Then 

< ^kA >= j a l ap{ HF}eM-F). (3.2) 
The first order contribution to this expression is 



Yl i°p4i J 2 P3 a P4exp(5^-/fc4ofc) 



P1P2P3P4. 



fpi ~l~ fp2 fp3 fpA 

u Pl + u P2 — u) P3 — uj P4 — ie 



(3.3) 



^PlP2P3P4 e ~ 0- 



It is equal to zero even with e > as the contractions will enforce f Pl + f P2 — f P3 — f P4 = 0. 
The second order contribution to ( |3.2|) is given by 

/ o,pa p ( — 2eOp 1 p 2 p 3 p 4 p 5 p a a pi a p2 a P3 a P4 a P5 a P6 + 

^ e ^piP2P3P4^P5P6P7Ps a p 1 a p2°'Pp3 a PA a 'p 5 a p 6 a 'P7 a P8) eX P( ^ ] _ fk a fc a k ) • 

k 

The integral involving O will be of the order of 1/e as the contractions will ensure the 
denominator in ( |2.4j ) to be equal to — 2ie. After rather long and tedious computations 
(which actually become quite simple with some practice) we obtain 

2e J a p a p fla^ aaaexp(—'^^ fka\ak) = Q p + P p — K p (3.5) 



' I am indebted to I. Kogan for pointing out this correspondence. 
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where 



P P— T7 7 f ~^^PlP2P3pJ( U Pl + U P2 - U P3 ~ W P 4 ) 

P1P2P3P4 JpJPiJPiJPaJPi 



(3.6) 



(fpi "I" fp2 fp3 fpi) [fippi "I" ^PP2 "I" ^PP3 "I" ^ 



PP4h 



n = V 4% A 2 

^P / ; f f f f f ^PlP2P3P4 

P1P2P3P4' /P ' /P1 ' /P2 ' /P3 ' /P4 



(3.7) 



5(o; Pl + w P2 - w P3 - uj p4 )(f Pl + f P2 - f P3 - f P4 Y 



and 



K P= ff f f X2 pP2P3pAfp + h~ f^ S ( UJ P +UJ P2 ~^P3 -Up A ). (3.8) 

P2P3P4 JPJPIJPSJP* 

The term Q p in ( |3.5|) came from contracting of a p with a p while the rest of the terms 
came from intercontracting those operators with the interaction. We recognize the kinetic 
equation ( |2.11|) in K p = while the term P p has the meaning of the entropy of waves^. 



As for the integral involving A 2 , its divergent part is determined by the formula 

1 TV 

S(x), e^O 



(x + ie)(x — it) e 
to give 

e Ja p a p Aa^a^aaAa^a^aaexp(—fa^a) = Q p + Pp. (3.9) 

The difference between (|3.5|) and (|3.9| ) gives us exactly the kinetic equation ( |2.11| ) as we 
expected. 

However the calculations we just completed imply more than just checking the validity 
of integrating by parts. We can use the above formulae to calculate the corrections to the 
wave density < a p a p >. We obtain 

< a p a p >= — + j a^a p (—£la)a^a)aaa-\ a) a) aaa) a) aa) exp(— / 'a) a) . (3.10) 

' Jp J 2 

These integrals have already been computed. The only difference with the integrals in 
( pT5| ) and ( pOl ) is the absence of the "vacuum diagrams" or Q p , and of course e will not 
be canceled. That makes a crucial difference with ( |3.4|) . Not only the expressions we shall 



t See ref. Q for details on the entropy definition. 
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obtain will be divergent with e — > 0, but also we need to extract the final parts from them, 
for example by using 

-5(x)--— + o(e) (3.11) 



(x + ie)(x — ie) e x dx 

We shall discuss the formulae like ( |3.11| ) later in more details. Now we just note that in 
addition to the divergent with e — > part, there will be a finite part which we shall not 
calculate here, leaving it to the discussion of the self energy corrections later in this paper. 
The answer will be 

< ala p >= - ^-(P p - K p ) + ^-P p + -0(e) = ^ + ^-K p + finite part (3.12) 

We see that the second order contribution to the wave density is clearly divergent and 
it should not be strange. The e regularization actually plays the role of cutting off the 
times greater than 1/e and if something is equal to 1/e, that means it is divergent with 
time going to infinity which can happen to the wave density if the kinetic equation is not 
satisfied. However if it is satisfied, we would rather expect no divergencies in ( |3.12| ). So it 
is reasonable to believe that the divergent higher order corrections to ( [3.12|) will sum up 
to be 1/e multiplied by the exact kinetic equation (|2.9|) . As long as the kinetic equation is 
satisfied, we can neglect those corrections. 



4. Lowest Order Corrections to the Kinetic Equation 

With the experience of calculating different correlation functions we acquired in the 
previous section, we are ready to start computing the corrections to the kinetic equation. 
All we need to do is to take ( |2.9|) and compute the imaginary part of the four point 
correlation function order by order by expanding exp(— F). The calculations are messy, 
but pretty straightforward. In this section we are are going to discuss some of the lowest 
order corrections to the equation ( 2.11 ), the ones which are not the corrections to the 
propagators, just to get acquainted with them. They consist of six terms. We begin with 
considering the first three of them. They are 

fpi fp2 ~ fp3 ~ JPi 



^ ] 87rAp lP2 p 3 p 4 Ap 3P4 p 5 p 6 Ap 5P6Pl p 2 

P1P2P3P4.P5P6 



fpi fp2 fp3 fp4 



■ X 



5{u Pl + u P2 - u P3 - u P4 ) f l 5+fp ° -+ (4.1) 

JpsJpeK^Ps u pe ~ U P3 ~ U P4> 
+ (3, 4) <- (5, 6) + (1, 2) <- (5, 6)) [d PPl + 

S pP2 $PP3 "PP4 ] 

(1,2) <-> (5,6) means that this expression differs from the explicitly written one by ex- 
changing the indeces. 
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Fig. 1: Some of the lowest order corrections to the four point function 

It is easy to represent the way how A's are connected with each other by a simple 
diagram borrowed from the (p A theory. The only feature we have to add is the little arrows 
on each line representing whether the line goes into the vertex or from the vertex. Since 
\1p2p3p4 i s included in any correction (it comes from the exact equation ( |2.9| ) rather than 
from the perturbative expansion), there is no need to represent it. The two remainings A's 
are represented by the first of the diagrams on fig. 1. 

As for the coefficient in front of the A's, it doesn't resemble anything in if 4 theory. 
The only familiar feature it retained from that theory is the product of all the propagators 
^-/(fpifp2fp2,fp4fpsfp&)- All the rest came from the interaction coefficients Q2.4] ) and ( |2~3| ) 
which are very nontrivial. The next section of this paper will be devoted to their analysis. 
Now we shall only note the fact that the first term in ( ^.1| ) coincides with the kinetic 
equation ( |2.11| ) itself multiplied by 



fps ^ fpe 



fpsfpe ( W P5 + U P6 ~ U P3 ~~ ^Pa) 



(4.2) 



The two intermediate lines in the diagram should be represented by (|4.2| ) . The second and 
the third terms in (|4.1|) are constructed from the first one, so we don't need to discuss 



them separately. 

The expression (|4.2|) has a simple physical meaning. The scattering represented by 
the diagram we consider can happen in two ways. Either the waves 1 and 2 scatter to 5 
and 6, and then 5 and 6 scatter to 3 and 4, or the waves 5 and 6 first scatter to 3 and 
4, and only then 1 and 2 scatter to 5 and 6. The first amplitude should be proportional 
to (715 + l)(riQ + 1) where n = 1/f is a wave density while the second to tl^uq. The 
"energy denominators", represented by the u>'s which are of the same origin as the similar 
denominators in quantum mechanical perturbation theory, will have opposite signs in these 
two cases so the amplitude will be proportional to (77,5 + 1)(^6 + 1) — ^5^6 ~ ^5 + ^6, that 
is exactly (gg). 
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The three terms corresponding to the second diagram from fig. 1 can be computed to 

be 

Eo 97r \ X \ I III + f p2 ~ f pa ~ ^ P4 v 

o^/i /\p 1 p 2 p 3 p 4 /\p 4 p 5 p 2 p 6 /\p 6 p 3 p ip5 I /\ 

V JP\JP2JP3JPA 



P2P3P4P5P6 



5{u Pl +u P2 -u V3 -u PA ) ^ fp ° -+ (4-3) 

JPS JP6 l W P5 + U P4 ~ U P6 ~ UJ P2 ) 

+(3, 1) <- (5, 6) + (2, 4) <- (5, 6)) [8 PPl + 

3pp2 "PP3 PP4 ] 

Here the intermediate waves propagate in the opposite directions, so the amplitude is 
proportional to (77,5 + l)ne — n^in^ + 1) ~ tiq — 77,5. 



5. Diagrammatic Technique 

We can proceed with calculating higher order diagrams in the same fashion. But first 
we shall consider another approach of calculating the corrections to the four point correla- 
tion function, the one which doesn't require the knowledge of the probability distribution. 
While it gives the same results, it is sometimes simpler to use for practical computations. 
This approach is based on the idea that the correlation functions should not depend on 
time. For example, if we compute the time derivative of the four point correlation function, 
we obtain 

< {H, a^a^a^a^} >= i Ka^a^a^a^ > [u\ + 0^2 — ^3 — ^4) + 

P1P2P3P4 

which expresses the four point correlation function in terms of the six point one. Clearly, 
the zero order value of the six point correlation function will give us the first order value 
to the four point one and so on. Quite analogously, to obtain the first order correction 
to the six point correlation function, we write that its derivative is zero to express it in 
terms of eight point correlation function and so on. In order to avoid the poles while using 
( |OD we replace ^2 u by u + * e - The sign of the e can be checked to be + by computing 
everything at finite time^. 

' That means by solving the equation for the time derivative of the four point function and 
choosing the imaginary part for the u) to be such that there's a finite limit at t — > 00. 
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Our aim is to find the perturbation series for the four point correlation function, as 
clear from the exact kinetic equation (|2.9|) . From the above consideration we see that the 
n-th order correction will be of the form 



, „ {-H"j. a j a 2 a 3 a 4} , 

Iff. Z^ ( " )+I£ \ 

t t 5>) + (n-l)ie ' 

i < a[a 2 a 3 a 4 > n =< ^7 \ , : >o (5.2) 

+ me 

where Hi is the interaction part of the hamiltonian (|1.1|) , 

P1P2P3P4 

and we retain just the zero order value for the correlation function on the right hand side 
of©. 



Let us see how ( |5.2| ) works. First we compute the Poisson bracket 



{Hi, a 1 a 2 a 3 a4} — i Xp^p^a^a^ap^a^a^a^ + . . . 

We represent the product o^a^a^a^ by a vertex with the incoming lines 1 and 2 and the 
outgoing lines 3 and 4. Taking the Poisson bracket means attaching another vertex to the 
original one with the lines pi, P2, P3, and 1, the last line going to the first vertex. Then, 



according to (|5.2|) we divide all we obtained by u\ + u>2 — W3 — 0J4 + it. Graphically, it is 
just the sum and difference of a/s corresponding to the incoming and the outgoing lines of 
the first vertex. Then we continue the procedure. Taking the Poisson bracket again means 
attaching another vertex to the diagram we are building, and then we need to divide by the 
sum and difference of w's corresponding to the first and the second vertex taken together. 
At last, when we arrive at the nth order picture, which is called the tree diagram, we need 
to compute the Gaussian correlation function, that means to close each outgoing external 
line of the tree with the incoming line and multiplying the expression by the product of 
the propagators 1/ / for each Gaussian contraction. Simultaneously our picture transforms 
itself into something resembling the standard ip theory diagram. The distinctive feature 
is the arrows on the lines (to distinguish between the incoming and outgoing lines) as we 
already saw in the previous section. Also, if we do everything honestly, we need to close 
the external lines 1, 2, 3, and 4 into one additional vertex, but that makes the drawing 
rather clumsy, so we shall be using the standard diagrams like those shown on fig. 1. 
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In short, we see that the right hand side of ( |5.2| ) is constructed by connecting A's 
from different Hi together by means of the Poisson brackets to form a tree diagram and 
multiplying u; in the denominator along the way, then by pairing all a with to form 
a real diagram writing the propagators 1/ / for each line which was formed by pairing and 
not by a Poisson bracket connections. The symbol u means first we take the sum of 
±o;'s belonging to the incoming or outgoing lines of the first vertex, then we take the sum 
of ±uj belonging to the incoming or outgoing lines of the first and the second vertex along 
the tree "trunk" combined and so on. 

We can reverse the procedure to say that first we should draw a <^> 4 diagram, then cut 
some of its lines to turn it into a tree diagram and write down the product of uj taken 
along the tree trunk. Then we multiply what we got by the product of / taken alone the 
trunk and divide it by the product of all /'s. And the last step should be summing over all 
different trees which can be curled up to form the same diagram. Each tree should have 
a direction along which we go and each trunk line then has a sign depending on whether 
this line has an arrow alone or opposite to that general direction. The sign of the diagram 
is the product of the signs of the trunk lines. It follows from the property of the Poisson 
bracket to have opposite signs in case when a is connected to or vice versa. 

We can go a little bit further in the formulation of our rules of reading diagrams. 
Namely if the tree have many branches, there are many ways we can sum up the w's. We 
can go along the first branch of the tree, then the second, or first along the second, then 
the first, or we can go up half of the first branch, then along the second, and then finish off 
the first one. And we need to take the sum of all the expressions involving u/s we obtain 
in this way. 

To compute this sum is a difficult job, but fortunately the expressions conspire in such 
a way as to give rise to the formula which can be written down using the Feynman rules 
given below which work for both the trees with branches and without, and so they give 
the general way of computing any diagram without doing the actual computation. 

The Feynman rules. 

In order to write down the expression corresponding to a given diagram, we need: 

1) To write down the product of A's corresponding to the vertices of the diagram, just like 
in </? 4 theory. 

2) To write down the product of all the propagators 1/ / corresponding to each line of the 
diagram, again just like in p 4 theory. 
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3) And to compute a special prefactor. First we need to make up a list of all possible 
expressions taken around each vertex of the diagram and around each combination 
of vertices. Then we need to cut some of the lines of the diagram to convert it into a tree 
diagram. The tree diagram by definition always starts from the external line while the 
external lines are thought as being connected together by some "external" vertex. The 
expression corresponding to the tree will be the product of all /'s taken along the trunk 
and the product of as many taken from the list we made up as there are /'s, 
each of the uj having one and only one particular omega belonging to one of the trunk 
lines. One can show there is one and only one subset of Y u with this property in the list 
we prepared. After doing that, we should replace all by XX^) ~~ ^ ne where n is the 
number of vertices we combined together to write down its J^w. n appears in front of e 
since in combining together ^(o;) — ie from adjacent vertices we have to sum their e as 
weil.t The last step here is summing up over all possible trees keeping in mind that the 
trees are taken with plus or minus sign depending on their intrinsic signs as was explained 
above. 

4) The analytical expression corresponding to a given diagram is the product of what we 
obtained in steps 1), 2), and 3). 

We would like to remark here that even though we obtained our rules by the method 
given in the beginning of this section, one can check in every particular case that the 
expansion of the probability density exp(— F) will give the same result. It is clearly so as 
one can show the equations like (|5.1j ) are always satisfied if the correlation functions are 
computed with the help of our probability density. Even though the probability distribu- 
tion possesses not only the four point vertices but also six and higher order vertices, they 
in fact conspire to give only the tp A theory diagrams as the higher order vertices are the 
products of the four point ones as is clear from (p73|), ( |2.4j) , etc. If one tries to rigorously 
prove the Feynman rules by expanding the probability density, one obtains the expressions 
which, though much more complicated than those given above, reduce to those simpler 
expressions after some tedious algebraic computations. It actually looks like a miracle that 
so many terms cancel in this procedure. To prove the algebraic identities appearing in this 

' We write here —it instead of it because we formulate these rules as computed in ( [2 .2D 
representation. The reader may easily check it by comparing the Y w given in the rules with the 
previous Y 00 which are computed together with the external vertex. 
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way is probably a very interesting but rather difficult task which we have not accomplished 
but which is not needed for our purposes. 

Let us compute the expression for the first diagram from the fig. 1 using these rules 
to see how they work in real life. The first tree we choose goes along the lines (1, 5). The 
expression it corresponds to is 

/1/5 



(uji -\- 002 — — 004 — 2ie)(uj5 + ujq — — 004 — ie) 

We have used a shorthand notation here f n = f Pn . 

Then there is another tree going along the lines (1,6). Its expression differs from 
only by replacing f 5 by fa. So we can sum them up together to obtain 

fi(h + h) 



(5.3) 



(5.4) 



(uji + U2 — ^3 — U4 — 2ie)(uj5 -f- ujq — <-o>3 — — ie) 

This is a part of a general rule. If there's a minimal loop consisting of the lines n and m 
going between two vertices, then we always have the factor of the form 

fn ~\~ fm f r r-\ 

(5-5) 

if the arrows of n and m are in the same direction and 

fn fm fr r\ 

(5-6) 

if they are in the opposite direction plus some other contributions which can appear if the 
tree trunk can go through both ends of the loop we consider without going along the lines 
n or m. 

Then there are trees going through (2, 5) or (2, 6). Their expressions are also almost 
the same as ( |5.3| ) with the replacement fi <-> f%. We sum them up to obtain 

(A + /2X/5 + /6) (57) 

+ UJ2 — ^3 — ^4 — 2ze)(cu5 -\- ujq — CU3 — CU4 — ie) 

Other trees without branches are (3, 5), (3, 6), (4, 5), and (4, 6). Their sum gives 

(/3 + /4X/5 + /6) (58) 

(uji + uj 2 - W3 - ^4 - 2ze)(wi + ix> 2 - u; 5 - uj 6 - ie) ' 

The remaining trees are with branches. These are (1, 3), (1, 4), (2, 3), (2, 4). Their 
expression can also be found using the Feynman rules. For example, (1,3) gives 

/1/3 



{uj x + u 2 — ~ ^6 — it) (^5 + uj 6 — u 3 — — ie) 
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(5.9) 



while all four trees can be summed up to give 

(/1 + /2X/3 + /4) 



(uji + 002 — — uj6 — ie) (^5 + ujq — UJ3 — — ie) 



(5.10) 



The minus sign appeared here as one of the lines of the tree trunks here goes the opposite 
way to the arrow of that line. 

The total expression for the diagram is the sum of (|5.7| ), ( |5.8| ), and ( |5.1U| ). If we want 



to compute the correction to the kinetic equation, we need to take its imaginary part. We 
obtain three terms proportional to the three possible delta functions. For example, the 
term with 5(u>3 + 004 — u)\ — U2) comes from (|5.7| ) and ( |5.8| ) and is proportional to 

(/1 + /2X/5 + /6) , (f3 + U)(h + fe) \ xf M * 

1 )0(UJ3 + UJ4 — UJi — OO2) = 

^3 + -UJq UJ 5 + UJ G - UJi - U 2 / (BU) 

Ui + f2-h-h){h + k) xf , , [ ' ' 

= 0(UJ3 + — UJ\ — UJ2) 

UJ3 + UJ4 — UJ5 — UJQ 



It coincides with the first term of Q4.1J ) after multiplying it by the necessary A's and 
by all the propagators 1/ /. 

There is an interesting observation applicable for the expressions we obtain after taking 
the imaginary part. The 8 functions which appear as soon as we take the imaginary part 
are always multiplied by the corresponding sum of /, that is all of them appear in the 
form S(J2 U ) J2f- It can be proved by analyzing the rules we discussed or, more easily, 
by noting that / = u should always be an exact solution of the kinetic equation ( |2l3| ) and 
S(^2u>) Yl 00 = enforces that last statement. 

As we start to compute higher order diagrams, additional feature we shall discover 
after taking the imaginary part is the appearance of the products of the odd number of 
5-functions in them. 

A very important property of the diagrams of almost any field theory which our 
diagrams are missing is the "block summing" property. Our rules are such that adding 
a block inside one diagram is totally different from having this block inside another one. 
This poses a major difficulty in dealing with our technique. 
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Fig. 2: Corrections to the propagators 



Now, we immediately ran into new difficulties as soon as we consider self energy 
corrections diagrams like those shown on fig. 2. Such diagrams have at least two lines with 
the same momentum (because of the momentum conservation at vertices) and, by using 
our rules, we arrive at divergent expressions of the form 

1 1 



uj — uj + ie ie 

(For example, the first diagram of fig. 2 has a vertex with a tadpole loop growing from it 
and its ^2(uj) is cleary zero.) It is in fact in the expressions of this kind that the coefficients 
in front of e become important. 

The appearance of the divergent expressions is nothing new; we already dealt with 
them in the section 2 to discover that they are actually zero by virtue of the kinetic 
equation. But there are still finite parts left. To compute them, we have to make use of 
the following important formulae: 

dx Ji± = ±^l[ dx ^M (5.12, 

fc /M '/j,/"-'W (5.13) 
[x ± ie) n n\ J x ± ie 

These formulae are correct if the integrals are convergent. If they are not, we have to 
include the boundary terms on the right hand side. 

Let us use them to calculate the first of two diagrams of fig. 2. We shall discover that 
this diagram, after adding to it all its generalizations, including those with many tadpole 
loops on the external lines, give the correction to the frequency uj in the kinetic equation 
( |2.11|) . Really, the first diagram of fig. 2 can be computed to give 

\2 fpi + fp2 ~ fpa ~ fP± X>(, , _|_,,, ,,, ,, \KlP5PlP5_ (K1A\ 

A P1P2P3P4 f f f f ^Pl + UJ V2 ~ ^P3 ~ UJ P4> f l 5 - 14 ; 

J Pi JP2 JP3 J PA JP5 
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8' appears as we have to use ( |5.12|) and ( |5.13|) in giving sense to the expressions we obtained 
using the rules discussed above while the convergence of the integrals in (|2.11| ) allows us 
not to worry about the boundary terms. 

The expression here, when combined with the kinetic equation (|2.11] ), clearly gives 
the correction to the frequency in the equation ( |2.11|) 

u p -> u p + ! ^-dq (5.15) 

While ( |5.14| ) is just the first term of the expansion due to ( |5.15| ) we can easily prove that 
( |5.15| ) can be obtained up to any order if one sums up the diagrams with all possible 
tadpole graphs on the external lines. However, more complex diagrams, like the second 
one from the fig. 2, can lead to the corrections which cannot be interpreted that easily. 

We would like to conclude this section with saying that the technique described here 
is more general than the standard field theory approach of the (p 4 theory. In fact, if we 
take the limit of f p — > u p we shall discover that the prefactor computed according to our 
rules goes into 1 and the whole expression turns into the standard </? 4 theory diagram. 
This should of course be expected as in this limit F — > H . So our Feynman rules should 
be treated as the "turbulent" generalization of the standard field theory rules. We believe 
many beautiful results are hidden in this new technique. 



6. Epsilon Expansion 



The kinetic equation ( [2.1 1|) is in general very difficult to solve (that is, to find such f p 



that it is satisfied) . However it has long been realized that if both A and uj are homogeneous 
functions of momenta, then we can solve that equation exactly. Namely, following ref. [|I| 
we choose 

UJ p =p° l , X P1P2P3P4 = X (PlP2P3P4)^U(pi,p 2 ,P3,P4)S(Pl +P2-P3 ~ Pa) (6.1) 

where U is a function depending only on the ratio of lengths of the momenta and the 
angles between them and Ao is a small constant. The parameter a is called the energy 
spectrum dimension, while f3 is the interaction dimension. Then the kinetic equation can 
be solved with the aid of the so-called Zakharov transformations, to give 

2 rv 

U = P\1= -[3 + dor 7 = -f3 + d~- (6.2) 
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d being the number of space dimensions. The latter of the solution corresponds to the flux 
of the wave number while the former to the energy flux. All that has been known for a 
long time. 

However the solution for f p we get in this way is just a first approximation to the full 
answer which should in fact be a series in powers of Ao obtained as a solution to the exact 
kinetic equation, with all the higher order corrections. The situation is reminiscent to 
that of the second order phase transition theory (on the theory of the second order phase 
transitions, see ref. || and references therein). The behavior of the perturbation series 
depends crucially on whether the physical dimension (measured in powers of momentum) 
of Ao is positive or negative. Let us denote that dimension by —k. Then, just from the 
analysis of the dimensions, we see that the series for f p should look like 

oo 

f P = J2 c - X oP 1+nK ( 6 - 3 ) 

c n being just some dimensionless coefficients. This is of course true only if the momentum 
cutoff does not appear explicitly in the series, that is all the integrals in the kinetic equation 
are convergent, but let us first examine this case. In general we are interested in the 
behavior of the spectrum f p at the distances much larger than the dissipation length. 
That means that we should actually take p — > when analyzing (|6.3| ) (like in the phase 
transition theory again). Now if k > then we can just retain the first term in the series 
( |6.3|) , while if k < 0, we can no longer do that. We must sum up the whole series and 
analytically continue it to the region where p — > 0. So, k plays the role of the e = d — 4 
parameter of the theory of the second order phase transitions. 



The method of phase transition theory to sum up the (|6.3|) -like series for k a small 
negative number is called the "e- expansion" (refs. and ||). We are going to apply this 
method to the theory we have here. 

First let us find what k is equal to. It is clear the physical dimension of the free wave 
part of the Hamiltonian ( |1 . 1| ) and the interaction part should be the same which leads us 
to 

a-<y = -K + f3-d + 2(-d - 7) + Ad 

or 

K = p + d-<y-a (6.4) 
We took into account that the dimension of the fields a and is —7/2 — d/2. 
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Let us note that upon substitution of (|6.2| ) into we obtain 

(3 (3 2 

k = a or k — a 

3 3 3 

which shows that actually k does not depend on the space dimensionality at all. This is 
different from what we usually have in the second order phase transition theory. 

Each next correction to the kinetic equation has a higher power of Ao- But the total 
physical dimension of all the terms should be the same. That means the total power of the 
momenta in each next term will be greater by k than that in the previous one. k > may 
lead to the divergent integrals in the perturbation series. However the divergency of the 
theory we consider does not depend on the power of the integrand only. It also depends 
on the properties of the interaction coefficient U. If U decreases exponentially when the 
ratio of the momenta becomes high, then all the integrals will be convergent. In this case 
the solution f p = p 1 is not modified at p — > for all k > (and even exact at n = 0) while 
K < is very difficult, if not impossible, to analyze. That is why, as was promised in the 
Introduction, we are going to consider the opposite case, U = 1, when the properties of 
the interaction do not influence the convergence of the integrals at afld. 

Now k > does lead to the divergent integrals. However it should not bother us as 
all the integrals are supposed to be cut off at some high momentum A where the energy 
dissipation takes place. The value of these integrals will depend on the value of the cutoff 
only and they change the kinetic equation ( |2.11 ) only by multiplying it by a number 



without changing its form. The situation is different when k < 0. Then the integrals in 
the perturbation series are more and more convergent and we should compute the sum 



There is also the third possibility, namely k = 0. In this case we can expect logarith- 
mically divergent integrals and the series (|6.3|) may turn into 



/ P =P 7 E c - A o lo S n (-) ( 6 - 5 ) 

n=0 ^ P ' 

The standard approach of the e-expansion theory tells us that if we are able to separate 
the most divergent diagrams which contribution are of the order of Aq log n (A/p) and sum 



' We must admit that it is the intermediate case which is usually fulfilled in Nature, that is U 
behaves as a power asymptotically, but we shall consider U = 1 to simplify everything as much 
as we can. 
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K(px,p 2 ,P^P4) = t t t t A p lP2 p3P 4 (/pi + /pa - /p 3 - /p4. 



them up, we can obtain a true asymptotic behavior f p for large logarithms. Now if k < 
we can apply the trick of replacing the logarithms in ( |6.5| ) by p K which will give us a 
modified wave density spectrum. This is a program we are going to fulfill. 
As was promised many times, we choose U = 1 in the interaction, or 

A PiP2P3P4 = A (P1P2P3P4)*<5(P1 +P2 ~P3 ~Pa)- (6.6) 

Also there will be an additional condition a < 1 we will have to impose later for the 
purpose of simplifying the calculations and which is very natural for the theory described 
by 0) (see §). 

Let us begin with analyzing the simplest bubble diagrams shown on fig. 1. The special 
form of the interaction (|6.6|) enables us to disentangle the momenta in it to bring it to the 
form 

~\~ T)^ (3 

P5P6 \P5 Pa P3 Pa ) 
where we denoted the kernel of kinetic equation by K, 

4n 

f f ~F f ' 'P1P2P3P4 \-> VI ' J V2 J P3 J P4/ 

JpiJp2JP3JP4 (6J 

S (u Pl + uj P2 — 0J P3 — oj P4 ) 

and assumed the integrals over p%, p2, ps and P4 will be computed later. 

The integral in (|6.7|) will be exactly logarithmically divergent at large momenta for 
K = 0. Namely if p§ and pe are very large we can safely assume pe = —p§ eliminating 
the 5-function in this way and be left with just one integral giving the log(A) at its upper 
limit. What is in the lower limit is not important, it is enough for us to write 

^ \u(p3,Pi) 

where u is some homogeneous function of the first order. The correction to the kinetic 
equation ( |2.11p we thus obtained is just 

,P2,2?3,P4)A slog 7^ ))(S PPl + 3pp2 fipP'3 ®PP4 ) (6-9) 

P1P2P3P4 \U{P3,P4)J 

Combining the correction we just obtained with the kinetic equation itself we arrive at the 
modified kinetic equation 

V K(p l ,p 2 ,p 3 ,P4)(l -2A slog ( A ^ ))($p Pl +S PP2 - S PP3 -S PP4 ) (6.10) 

P1P2P3P4 N ' 
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s being a number coming from the angular integration. 

It is not difficult to persuade ourselves that the sum of the contributions of the second 
and the third terms of (|4.1|) will be equal to (|6.9|) . Really, it is evident that 

^ ] 87rAp lP2 p 5 p 6 Ap 3 p 4 p 5 p 6 Ap 3P4Pl p 2 x 



P1P2P5P6P3P4 

[5pp 3 + 5pp 4 ] 



(fpi + fp 2 - fp 5 ~ fp e )^( U Pl + U P2 ~~ U P5 ~ U p 6 )(fp3 + fp4 



fpifp2fp3fp4,fp5fp6( U! P3 JrU, P4 ^Ph W P6 ) 

since the expression above is antisymmetric under the interchange of pi,p2 with P5,Pq. 
On the other hand, the same expression with 5 PPl + 5 PP2 will give us one half of (|6\9|). A 
similar thing can be said about the third term of ( |4.1| ) . Together all of them give us twice 
the (|Op. 

Unfortunately, these are not the only divergencies associated with this diagram. It 
is clear the diagram is going to be divergent for p*, — > or pg — > 0. The degree of the 
divergence is just (3/2 — 7, so it is —a/2 for the energy spectrum when, due to the condition 
K = 0, /3 should be equal to 3a or logarithmic for the wave number spectrum when (3 — 2a. 
This divergence will bring additional correction to the kinetic equation 

q 2 

— or \og(q/a) (6.11) 
a 2 

where a is the IR cutoff while q is the momentum of the order of the external momenta. 

However, if q is sufficiently small, we can neglect the IR term we just obtained in 
comparison with the UV one. We admit that this statement is rather vague. In fact, to 
have a satisfying proof, we need to sum up all the IR corrections corresponding to different 
diagrams. Most probably, there is some sort of condensate being formed like in the very 
well known case of Bose gas which will cancel the divergencies. But this has not been 
shown yet. We intend to devote another paper to the IR problem. Here we just assume we 
can neglect all the IR terms. This simplifies things a great deal as we can neglect the IR 



corrections to the frequency as well. The frequency correction (|5.15|) is IR divergent with 
the same power as the bubble diagram we considered. We shall neglect that correction 
as well while leaving the task of evaluating the physical effects it can lead to to another 
paper. Acting consistently, we also neglect all the self energy correction diagrams. 
Now we consider the second diagram given by 

dp 5 dp 6 7 7 . 5 6 ^(P5Pe)^ 8(p4 - p 2 + P5 -Pe) (6.12) 

P5JW5 -Pe +P4 -P2) 

It turns out the condition a < 1 makes it convergent. Really, p$ ~ pq in this case, thus 
enforcing p§ — p§ << 1. The power of the integrand is then estimated to be 7 — 1 — 27 + 
f3 + d=-l+a<0. 
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Fig. 3: The most important logarithmically divergent diagrams 



The nontrivial statement about more complex diagrams is that when we go to the 
higher orders, only the diagrams shown on fig. 3 which are the direct generalization of 
( f4.1|) will give us the necessary power of the logarithm. The expressions corresponding to 
them can be easily established by using the rules of the previous section to be just the 
product of ( |4.2|) for each loop. In addition to the logarithm to the power of n each of such 
bubble diagram will bring a factor of (— 2s) n (n + 1). n + 1 appears as the diagrams come 
with 2n replicas different by the original ones by the interchange of indices which give us 
the same expression as the original ones divided by 2. 

We can simply sum up the whole series 

oo / A \ 1 

Y(n + l)(-2s\ ) n log n ( -— — ) = ? (6.13) 

So we obtained the full kinetic equation, with the most important higher order con- 
tributions summed up 

3pp2 3pp'i 3 PP4 ) 

^PlP2P3pAfpi + fp2 ~ fps ~ fpJdj^Pl + U P2 ~ U P3 — U Pi) _ n (6-14) 

2 



1 + 2sA log ( M( A~ 4) ) f Pl f P2 f P3 f P4 



All of the remaining diagrams will not contribute to the renormalized kinetic equation 
because each of the more complex diagram will have the lines going in the opposite direction 
and we shall always encounter expressions like ( |6.12| ) while calculating them. 

Let us begin with the integrating out the momentum which goes around the minimal 
loop (and the minimal loop is the one which goes just between two vertices). Then we shall 
obtain the diagram where this loop is replaced by a dot multiplied by the logarithm of A. 
Really, the expression for the diagram came from the sum over different trees. If the tree 
went alone one of the propagators of the loop, then the the loop momentum should enter 
the expression in the form of a multiplier ( |5.5[ ) or ( |5.6| ) depending on whether this loop 
has parallel or opposite arrows. The multiplier ( |5.5|) is integrated out to give log (A) while 
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the ( |5.6|) gives the finite answer. What is left after the integration is exactly the very same 
diagram with that loop shrunk to a point. As for the trees not going through the loop, the 
integration over the loop momentum gives the convergent integral as those trees do not 
give rise to the inverse loop propagator / in the numerator (we remember that the /'s in 
the numerator are taken from the trunk of the tree) thus making the integral convergent. 
If we continue this procedure shrinking the minimal loops one by one, sooner or later we 
encounter the loops with the opposite arrows Q5.6| ) and the integral will not acquire the 
necessary power of the logarithm. The only exception to this case is the diagrams shown 
on fig. 3. That is why they are the only diagrams we need to sum up. Of course, there 
is an IR problem left. After integrating over the minimal loop, we are left with not only 
the logarithms, but also with the positive powers of the external momenta due to the IR 
terms (|6.11| ) (or additional logarithms if the system transports waves instead of energy). 
They can lead to divergent integrals in subsequent integrations. For now, we just neglect 
those IR terms. 

Now our task is to find such f p which solve the equation ( |6.14| ). We should look for 
f p in the form of (|6.5| ). If we plug it into (|6.14| ), we can discover that the values of the 



coefficients c n depend crucially on the form of the interaction (except for c±). Fortunately, 
it is not the exact form of f p that we are looking for. We are interested in the asymptotic 
behavior of f p for large log(A/p). This behavior can be shown to be universal. It is 
more or less clear even without any computations. ( |6.14| ) has an additional log 2 (A) in the 
denominator which should be canceled by the additional log in /. That implies that the 
asymptotic behavior of / is just 

/ p =^log-i (£\ (6.15) 

The proof of ( |6.15| ) is rather messy and not instructive at all, so we shall not give it 
while noting that it exists. 

We can summarize the result obtained as follows. If a < 1, the interaction is chosen 
in the form (45), and k = 0, then the asymptotic behavior of the inverse wave density f p 
is given by (|6.15|) . If k is a small negative number, then the asymptotic behavior of f p will 
be 

f P =P 1 ~^ (6.16) 
and the energy spectrum density will be 

E p = V a -~ i+ -* K (6.17) 
The spectrum was modified by Ik which is called the anomalous dimension. 
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7. Conclusions and open problems 

So we have shown how one can deal with the turbulence correlation functions for the 
system with a small interaction. There is a task of classifying the IR contributions to the 
self energy and other corrections not accomplished here, but it seems the technique we 
developed is powerful enough for this work to be done in the nearest future. 

However, the most challenging problem to solve would be to find the correlation 
functions for the systems described by the Hamiltonian 



P1P2P3P4 



With a particular choice of A this Hamiltonian describes the motion of the incompressible 
fluid (see ref. [JJ). This theory possesses the integrals F in the very same fashion as the 
theory ( |1 . 1| ) . Unfortunately we know of no way to find them explicitly as the equations 
of motion which follow from (|7.1|) are impossible to solve. Still some of their properties 
can probably be studied by extrapolating the properties of ( |2.2| ) even though the most 
pointblank extrapolation, that is uj — > 0, is most probably not possible as a phase transition 
may take place. 

It would also be especially interesting to study these probability distributions in two 
dimensions. We may expect them to give rise to certain conformal field theories with un- 



usual properties. The distribution (|2.2|) is essentially nonlocal even for the initially local 
Hamiltonian while the actions considered in conformal field theory are usually supposed 
to be local. And while the distribution (|2.2|) is explicitly real, the theory it describes can 
naturally be nonunitary because of the possible anomaly in the wave number conservation. 
The idea of the conformal invariance of the turbulence correlation functions was put for- 
ward and exploited in Q for the case of hydro dynamic turbulence but perhaps the study 
of the differences of the standard statistical distributions and the turbulent distributions 
can shed more light on this subject. It is also perhaps possible to formulate a simple 



statistical model like Ising model with probability distributions similar to (|2.2| ) and study 
its properties by direct computation. 



The author is grateful to A.M. Polyakov for sharing his ideas and especially for his 
support and interest in this work and to I. Kogan for many very useful discussions. 
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